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THE  DERIVATION  OF  SIMPLE  POLES 
IN  A TRANSFER  FUNCTION 
FROM  REAL  FREQUENCY  INFORMATION 


Abstract 


This  report  is  the  third  in  a 
series  of  three  that  evaluate  a 
technique  (frequency-domain  Prony) 
for  obtaining  the  poles  of  a transfer 
function.  The  main  objective  was  to 
assess  the  feasibility  of  classifying 
or  identifying  ship-like  targets  by 
using  pole  sets  derived  from  frequency- 


domain  data.  A predictor-correlator 
procedure  for  using  spectral  data  and 
library  pole  sets  for  this  purpose 
was  developed.  Also  studied  was  an 
iterative  method  for  reducing  the 
effects  of  noise  and  a technique 
based  upon  magnitude-only  spectral 
data.  ^ 


Introduction 


A technique  for  using  real- 

frequency  data  to  obtain  the  simple 

poles  of  a transfer  function,  and  the 

results  of  its  application  to  real 

electromagnetic  (EM)  data  were 

described,  respectively,  in  two  pre- 
1 2 

vious  reports.  ’ The  procedure  was 

referred  to  as  frequency-domain  Prony 

because  of  its  close  similarity  to  a 

method  of  extracting  poles  from  time- 

domain  waveforms  based  upon  Prony 's 
3 

technique.  First  applied  to  spec- 
ified  pole  sets  to  evaluate  its  char- 
acteristics'*" it  was  then  used  for 

real  EM  data  to  provide  a more  rel- 

2 

evant  test  of  its  potential.  Some 
variations  of  the  basic  technique 
were  also  studied  and  described.  In 


this  work,  our  main  objective  was  to 
assess  the  feasibility  of  classifying 
and/or  identifying  radar  targets  by 
use  of  pole  data  derived  from 
frequency-domain  information.  We 
assumed  the  problem  to  be  one  of 
estimating  whether  an  unknown  object 
belongs  to  a class  of  concern 
(classification)  and  if  so,  which  one 
of  that  class  it  is  (identification) . 
After  obtaining  some  pole  results 
derived  from  frequency-domain  Prony 
for  simple  ship-like  objects,  we 
attempted  to  use  such  pole  sets  to 
establish  which  target  from  a library 
of  pole  sets  is  most  like  an  unknown 
target  whose  scattered  field  spectrum 
has  been  obtained. 
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Results  pertaining  to  1)  a 
"shifting"  procedure  applied  to 
noisy  data;  2)  a vari-linear 
frequency-domain  predictor;  and 


3)  a pole  extraction  technique 
based  upon  magnitude-only  data 
are  discussed  in  the  text  and 
appendices . 


Pole-set  Results  of  Simple  Ship-Like  Models 


In  the  first  part  of  this  study, 

before  frequency-domain  Prony  had 

been  developed,  pole  sets  were  derived 

from  applying  time-domain  Prony  to  the 

computed  transient  waveform  of  several 

ship-like  targets  of  increasing  com- 
4 

plexity.  Because  of  various  uncer- 

4 

tainties  encountered  in  this  effort, 
it  was  deemed  necessary  to  consider 
less-complex,  targets  while  still 
retaining  some  features  characteristic 
of  ships. 

Consequently,  the  series  of  simpler 
targets  depicted  in  Pig.  1 were 

•k 

modeled  in  the  frequency  domain,  and 
their  backscattered  fields  were 
obtained.  Frequency  spectra  were 
calculated  for  two  angles  of  inci- 
dence, 30°  and  60°,  with  respect  to 
the  long  axis  of  the  wire  models 
(Fig.  2).  The  straight  wire  was 
included  for  reference  purposes  at  an 
angle  of  60°  only.  The  results  of 
applying  frequency-domain  Prony  to 
these  spectra  are  summarized  graphi- 
cally and  numerically  in  Fig.  3. 


Since  frequency-domain  Prony 
had  been  developed  in  the  interim. 


We  observe  that:  1)  the  greatest 

variations  in  pole  values  obtained 
for  the  two  angles  of  incidence  occur 
in  their  real  (damping)  components; 

2)  quite  good  agreement  exists  in  the 
imaginary  components  of  most  of  the 
two  pole  sets  for  the  same  target; 

3)  some  poles  do  not  appear  in  both 
sets,  even  though  their  invalidity  is 
not  indicated  by  either  a small  resi- 
due value  or  large  negative  real  com- 
ponent; 4)  the  smooth  pole  trajectory 
of  the  straight  wire  is  changed  con- 
siderably by  addition  of  the  wires 
used  in  models  1-3;  and  5)  significant 
differences  are  also  discernable  in 
the  backscatter-f ield  spectra.  These 
results  demonstrate  that  the  pole  sets 
do  contain  potentially  useful  infor- 
mation for  classifying  and/or  iden- 
tifying radar  targets  but  not  without 
some  uncertainty  caused  by  possible 
discrepancies  arising  in  the  pole 
extraction . 

The  actual  poles  of  a given  object 
are  uncertain  primarily  for  two  rea- 
sons. First,  without  knowing  apriori 
what  number  of  poles  will  contribute 
to  the  transfer  function  over  a given 
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frequency  interval,  we  over-estimate 
the  number  in  the  data  when  using 
frequency-domain  Prony.  This  produces 
some  curve-fitting  poles  which  may  not 
all  be  recognizable.  Second,  some 
pole  locations  may  vary  due  to  numer- 
ical inaccuracy  as  the  problem  param- 
eters (e.g.,  angle  of  incidence)  are 
changed,  while  others  may  totally 
disappear  from  the  data  because  they 
are  not  excited  for  a given  set  of 
parameters.  Noise,  while  not  a factor 
in  the  computed  data  presented  here, 
can  also  be  responsible  for  such 
effects. 

One  approach  to  this  problem  that 
has  worked  especially  well  for  noisy 
data  in  the  time  domain,  should  be 
equally  useful  for  reducing  the  effects 
of  parameter  variations  as  well.  By 
varying  the  parameters  of  the  extrac- 
tion process  we  can  derive  several 
sets  of  poles  from  the  same  data  and 
then  plot  them  on  a single  graph. 

When  this  is  done,  the  valid  poles 
are  found  to  occur  in  clusters,  and 
their  positions  can  then  be  estimated 
by  computing  the  mean  coordinates  of 
each  cluster.  In  addition,  the 
standard  deviation  of  the  points  in 
each  cluster  can  help  to  indicate  the 


uncertainty  in  the  location  of  that 
particular  pole.  Although  not  devel- 
oped sufficiently  for  use  in  this 
study,  some  preliminary  results  of 
utilizing  this  technique  with 
frequency-domain  Prony  are  described 
in  Appendix  A.  Our  previous  obser- 
vation that  the  greatest  variance  in 
pole  location  usually  occurs  in  the 
damping  component  is  confirmed  by  the 
results  of  Appendix  A and  theoretical 
analysis  as  well."* 

It  is  interesting  to  compare  the 
lowest  frequency  poles  for  the  four 
objects  studied  above.  While  the 
straight  wire  has  a ju)-axis  component 
of  0.792  MHz,  a lower  value  of 
0.758  MHz  occurs  when  an  appendage  is 
added  near  one  end  of  the  wire  (ship- 
model  1).  However,  when  the  wire 
appendage  is  added  at  the  center  (as 
in  ship-model  2) , the  resonance 
(0.788  MHz)  is  only  slightly  lower 
than  the  wire-only  value  because  the 
appendage  has  little  influence  on  the 
even  current  distribution  character- 
izing the  lowest-f requency  resonance. 
For  ship-model  3,  the  (average)  reso- 
nance (0.751  MHz)  is  the  lowest  of 
the  four  because  of  the  combined 
effects  of  both  appendages. 


Use  of  Pole  Sets  for  Object  Classification 
and/or  Identification 

GENERAL  the  poles  of  radar  targets  is  their 

A major  motivation  for  studying  potential  use  for  classification 
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and/or  identification  (C/I)  independ- 
ent of  aspect  and  polarization.  As 
we  have  seen  above,  an  appearance  of 
aspect  dependency  can  occur  because 
some  pole  responses  may  not  be  excited 
at  certain  aspects  and  the  computed 
locations  of  those  that  are  excited 
may  vary  slightly  with  the  problem 
parameters.  This  indicates  an  inher- 
ent uncertainty  of  indeterminant  mag- 
nitude no  matter  how  the  poles  or 
information  derived  from  them  might 
be  used.  Noise,  of  course,  can  be 
expected  to  increase  this  uncertainty. 
In  any  case,  it  would  be  expected  that 
a common  element  of  any  pole-based  C/I 
approach  would  be  a pole-set  library. 

POLE-SET  CORRELATION 

Perhaps,  the  most  straightforward 
way  to  use  the  pole  sets  in  a C/I 
exercise  would  be  to  correlate 
directly  the  poles  obtained  from  the 
scattered  field  of  an  unknown  target 
with  each  of  the  pole  sets  in  the 
library.  This  might  be  done  opti- 
cally, for  example,  by  the  use  of 
transmission  filters  where  the 
target-pole  set  is  overlayed  with 
each  library  set,  and  the  relative 
correlation  is  established  in  propor- 
tion to  the  amount  of  light  trans- 
mitted. A digital  correlation  tech- 
nique could  also  be  easily  conceived, 
but  this  approach  presupposes  that  a 


pole  set  will  be  derived  from  the 
spectrum  of  the  unknown  target. 

Because  we  have  found  the  pole 
extraction  procedure  may  be  numeri- 
cally inconsistent,  computation  of  the 
pole  set  from  measured  data  is  not 
always  advisable.  In  addition,  using 
the  Prony  technique  requires  substan- 
tial computing  capability,  and  without 
access  to  a large  computer,  a near 
real-time  target  C/I  may  not  be 
achievable.  For  these  and  other  rea- 
sons, the  direct  use  of  the  poles  for 
C/I  may  not  be  the  most  useful 
approach  to  take. 

A PREDICTOR-CORRELATOR  APPROACH 

An  alternative  has  been  tested  in 
£ 

a transient  format,  which  involves 
the  concept  of  a linear  predictor  of 
the  form 

f(t)  = £ t). 

n=l 

t _>  N A t,  (la) 

where  N is  the  number  of  poles  in  the 
response.  The  predictor  provides  a 
way  to  obtain  sampled  values  of  a 
transient  signal  for  all  t > T in 
terms  of  the  given  samples  in  the 
interval  t < T,  and  a set  of  coeffi- 
cients a.,  which  are  defined  by 
N,n  1 
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a=l 


s At 
a s 


N 

s 


+ X 


a.,  s 
N,ot 


N-a 


(lb) 


where  is  the  ath  complex  resonance 
(or  pole) . Note  that  for  the  predic- 
tor to  be  applied  to  a given  f(t), 
the  set  of  poles  for  the  object  that 
produced  f(t)  must  be  available. 

Suppose,  however,  that  the  object 
is  unknown  and  to  be  determined.  If 
a library  of  pole  sets  is  available, 
where  one  may  correspond  to  the  object 
from  which  f(t)  has  been  obtained, 
then  the  coefficients  of  Eq.  (lb)  for 
each  library  pole  set  could  be  tried 
successively  in  Eq.  (la).  In  each 
case  the  known  signal  f^(t)  would  be 
used  on  the  right  hand  side  of  the 
equation,  which  would  lead  to  a set 
of  computed  signals  f (t) . Upon  com- 
parison of  fc(t)  with  f^(t)  to  estab- 
lish the  pole  set  that  minimizes  their 
difference,  the  reference  target  most 
like  the  unknown  could  be  identified. 
This  procedure  would  thus  allow  an 
identification  to  be  made  without  the 
necessity  of  finding  the  poles  from 

fk(t).6 

The  predictor  just  discussed  is 
based  on  transient  data  while  we  are 
interesed  in  the  use  of  spectral  or 
frequency-domain  information. 

Although  a transient  signal  could  be 
derived  by  Fourier  transforming  the 


spectral  data,  such  an  extra  step 
seems  unnecessary.  It  would  be  more 
direct  to  confine  the  process  entirely 
to  the  frequency  domain  and  find  a 
frequency-domain  equivalent  of  the 
linear  predictor. 

Although  a direct  frequency-domain 
analogue  of  Eq.  (1)  does  not  seem  to 
exist,  a procedure  can  be  identified 
that  operates  in  a strikingly  similar 
way.  To  appreciate  this,  let  us  first 
examine  the  application  of  the  time- 
domain  linear  predictor  [Eq.  (1)]  in 
more  detail.  First,  it  requires  a 
certain  minimum  amount  of  information 
before  it  can  be  implemented.  For  an 
N-pole  system,  the  N-pole  values  are 
needed,  which  amounts  to  2(N/2)  = 
N-real  numbers  since  the  poles  occur 
in  complex  conjugate  pairs.  Second, 
the  N-resonance  amplitudes,  or  their 
equivalent,  are  also  required  to 
specify  the  system  fully.  This  intro- 
duces additional  N-real  numbers  for 
a total  of  2N.  In  Eq.  (1),  these 
resonance  amplitudes  are  not  explic- 
itly involved  but  enter  implicitly 
instead  through  the  initial  samples 
of  f(t).  Thus,  from  these  2N-real 
numbers  we  can  obtain  all  other  values 
of  f(t). 

Now  let  us  consider  what  is 
required  to  develop  an  analogous  pro- 
cedure in  the  frequency-domain. 

Because  the  same  amount  of  informa- 
tion is  involved,  2N-real  numbers 
will  still  be  required.  Again,  this 
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might  include  the  N-pole  values,  given 
by  N-real  numbers,  while  the  other  N 
numbers  could  come  from  the 
N-resonance  amplitudes.  Alterna- 
tively, as  in  Eq.  (1),  they  could 
come  from  samples  of  F(to).  From  this 
data,  other  F(w)  values  could  be 
derived  for  any  w value  in  the  fre- 
quency range  over  which  the  2N-real 
numbers  were  obtained.  In  this  sense, 
a frequency-domain  linear  predictor 
differs  from  its  time-domain  counter- 
part, which  applies  for  all  subsequent 
time.  Of  course,  if  all  system  poles 
were  available  and  used,  then  the 
frequency-domain  version  would  also 
apply  for  all  frequencies.  Note  also 
that  successful  application  of  the 
time-domain  version  does  require  that 
all  poles  contributing  to  f(t)  be 
included  in  the  linear  predictor,  so 
that  use  of  the  linear  predictor  for 
all  time  also  requires  all  contrib- 
uting poles. 

A frequency-domain  linear  predictor 

could  be  written  in  the  same  way  as 

Eq . (1);  whereby,  all  F((jj)  for  lx)  > u)^ 

are  given  in  terms  of  F(co)  samples 
* 

for  a)  < 0),,.  However,  this  is  unnec- 
essarily  restrictive  since  we  can 
generally  predict  and  sample  at  rela- 
tively arbitrary  frequencies  over  the 
spectrum  of  concern.  Therefore,  we 
simply  use 
* 

This  is  shown  in  Appendix  B. 

^The  actual  ti)  used  must  satisfy 
the  appropriate  sampling  criteria. 


N 

F(jw)  = Ra/(jw  - sa).  (2a) 
a-1 

where  the  corresponding  time-domain 
expression  is  given  by 
N 

r(t)  - Z V • <2b) 

oi=l 

Proceeding  as  for  the  time-domain 
linear  predictor,  assume  that  the 
N-complex-conjugate  s^  are  available. 
Then  given  N-real  numbers  from  samples 
of  F(w),  Eq.  (2a)  provides  the  basis 
for  computing  the  as-yet  unknown  R . 
Suppose,  for  example,  we  are  given 
measured  real-frequency,  complex  sam- 
pies  of  F(jo)k)  = M^,  k = 1, . . . ,N. 
Then  from 

\-£  Xk1>Ea1>-  <3‘> 

- a 
a=l 

compute  the  vector  where 


.(i) 


(3b) 


any  N of  the  2N-data  points  provided 

by  the  N-complex  samples  M(jaj^)  may 

be  used,  and  the  superscript  i denotes 

the  library  pole  set  being  employed. 

* 

It  should  be  noted  that  R = R 

* * a a+N/2’ 

sa  = sa+N/2  and  \ = \+N/2*  If  this 
process  is  performed  for  each  of  the 

Frequency-domain  Prony  currently 
requires  both  phase  and  amplitude 
data  in  general;  see  Appendix  C. 
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M-pole  sets  in  the  library,  then  M 
vectors  R^  , i=l,...,  M will  be 
obtained.  They  by  using  the  remaining 
N data  points  provided  by  the  measured 
results  F(jco^),  compute 


where 

^ ■ I 'f’/O”  - 4l))-  <4b) 

a=l 

and  M^^  and  its  complex  conjugate 
are  both  used. 

We  refer  to  the  overall  process  as 
a predictor-correlator.  The  subscript 
xi  denotes  R as  a measure  of  simi- 
larity between  the  unknown  target  x 
and  library  entry  i.  The  smaller  R ^ 
is,  the  greater  the  similarity,  and 
vice  versa.  A normalized  version  of 

R . that  varies  between  0 and  1, 
xi 

where  the  latter  denotes  perfect  cor- 
relation, is  perhaps  more  conven- 
tional, and  we  define  it  by 


xi 


1 + R 


xi 


(4c) 


Thus,  by  using  part  of  the  measured 

data  to  find  a residue  vector  R^ 

a 

corresponding  to  each  pole  set,  we 
can  predict  the  response  at  the  other 
measured  points  and  compare  it  with 
the  actual  values.  The  results  giving 
best  agreement  with  the  measurements 
(smallest  R ) establish  that  the 


pole  set  from  which  they  were  obtained 
is  most  like  the  unknown  target,  with 
statistical  analysis  necessary  to 
establish  confidence  and  probability 
bounds.  The  advantages  of  this 
approach  are  several: 

• No  pole  sets  need  be  derived 
from  measured  data  obtained  in 
an  operational  mode. 

• The  library  pole  sets  required 
can  be  established  under  more 
favorable  laboratory  conditions. 

• The  excitation  of  all  poles  is 
unnecessary  for  the  predictor- 

correlator  to  work. 

• The  entire  operation  is  per- 
formed in  the  frequency  domain, 
and  generating  time-domain 
information,  either  via  trans- 
forms or  directly,  is  unnec- 
essary. 

• It  is  well  suited  for  analysis 
from  an  information-content 
viewpoint  thus  providing  an 
opportunity  for  better  under- 
standing and  more  efficiently 
utilizing  the  data. 

Some  numerical  results  obtained 
from  Eqs.  (3)  and  (4)  are  summarized 
in  the  next  section. 

SOME  NUMERICAL  RESULTS 


The  procedure  discussed  above,  as 
given  by  Eqs.  (3)  and  (4),  was  applied 
to  some  of  the  ship-model  data  previ- 
ously presented  in  Figs.  2 and  3.  For 
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this  test,  the  backscattered  field  of 
each  ship-model  was  used  successively 
as  the  unknown  target  return,  while 
the  library  consisted  of  the  pole 
sets  of  all  three  ship  models.  We 
used  poles  obtained  for  an  incidence 
angle  0 = 60,°  while  back  scattered 

data  for  0 = 60°  (case  A)  and  30° 

(case  B)  were  employed  for  the 
predictor-correlator.  The  results 
are  summarized  in  Table  1. 

Some  representative  graphical  data 
which  depict  the  steps  involved  in 
implementing  Eqs.  (3)  and  (4)  are 
shown  in  Fig.  4 for  Case  A.  Plotted 
there  is  the  total  imput  data  spectrum 
for  ship  model  1,  as  well  as  those 
points  used  for  the  residue  calcula- 
tion [Eq.  (3)]  and  the  predictor- 
correlator  [Eq.  (4)].  Note  that 


Table  1.  Target  correlation  matrices 

R .. 
xi 


Case  A 

x/i 

A 

R . 

XI 

1 

0D  t =60°; 

Poles 

2 

°RCS  " 6°° 
3 

1 

1.000 

0.329 

0.637 

2 

0.316 

1.000 

0.701 

3 

0.523 

0.536 

1.000 

Case  B 

x/i 

1 

0D  - =60°; 
Poles 

2 

eRCS  " 3°° 
3 

1 

0.760 

0.330 

0.389 

2 

0.235 

0.561 

0.  340 

3 

0.461 

0.531 

0.736 

these  two  sets  interleave  in  fre- 
quency, i.e.,  every  other  data  point 
is  used  for  the  residue  calculation, 
with  the  in-between  points  subse- 
quently employed  for  the  predictor- 
correlator.  Following  these  plots 
are  graphs  that  pertain  to  each  pole 
set  and  show  the  predicted  response 
and  the  difference  in  magnitude 
between  the  predicted  and  actual 
result,  |M(i)  - M I.  The  latter, 
specifically,  is  useful  in  showing 
where,  in  the  frequency  spectrum 
covered,  the  differences  between 
and  are  largest , and  thus  where 
the  discriminability  between  the 
various  targets  may  be  greatest.  The 
possibility  for  establishing  a measure 
of  information  content  per-unit- 
frequency  interval  with  respect  to 
maximizing  differences  in  target 
observables  along  these  lines  is 
interesting  and  deserves  further 
examination. 

Results  for  the  correlation  study 
are  summarized  in  Table  1.  We  observe 

A 

the  R. . has  a value  of  unity  for  all 
JJ 

three  cases  (i.e.,  j = 1,  2 or  3)  in 
which  the  same  0 is  used  for  the 
pole  and  field  calculation  (Case  A, 
Table  1) . This  shows  that  the  pre- 
dicted return  matches  that  measured, 
when  the  pole  set  which  corresponds 
to  the  target  is  used  in  Eq . (3).  We 
expected  this  because  the  same  target 
data  was  used  to  compute  the  respec- 
tive pole  sets  and  the  correlation 
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number  . . The  corresponding  R , 
i ^ j values  for  this  case,  however, 
exhibit  values  that  vary  between  ~0.3 
and  ~0.7,  which  establishes  a measure 
of  the  dissimilarity  between  the 
three  ship  models.  Note  that  the 
smallest  R, . values  occur  when  the 

ij 

single-masted  models  are  compared, 
which  is  a reasonable  result  when 
the  geometrical  differences  between 
the  three  models  are  considered. 
Correlation  of  the  two  one-mast 
models  with  the  two-mast  model  re- 
sults in  the  larger  R. . values. 

ij 

Again,  this  agrees  with  what  we 
might  expect. 

When  the  data  used  for  the  pole 
calculation  are  different  from  those 
to  which  the  predictor-correlator  is 
applied,  then  no  R values  are  unity 
(Case  B,  Table  1).  Evidently,  this 
is  because  of  the  slight  variation  in 
the  poles  which  characterizes  the 
different  incidence  angles,  and 
occurs  from  current  numerical  limita- 


tions (as  discussed  earlier) . How- 

A 

ever,  we  still  find  that  R. . for  each 

jj 

target  spectrum  is  largest.  Presum- 
ably, by  developing  an  internally 
consistent  averaged  pole  set  for  each 
library  entry,  such  angular-dependent 
variations  could  be  reduced. 

Although  the  technique  presented 
above  appears  to  have  some  potential 
for  target  C/I,  considerably  more 
effort  would  be  necessary  to  make  it 
practical.  The  effects  of  noisy  data, 
aspect  variations,  accuracy  of  the 
library  pole  sets,  etc.,  all  need  to 
be  considered.  In  addition,  other 
ways  of  computing  R . , for  example 
using  the  step-wise  predictor  de- 
scribed in  Appendix  B as  an  alterna- 
tive to  the  approach  described  here, 
should  be  evaluated.  More  funda- 
mentally, the  information  content  of 
the  available  data  should  be  explored 
as  a means  of  establishing  data  re- 
quirements, error  estimates,  C/I  con- 
fidence, etc. 


Summary  and  Conclusions 


This  report  is  the  third  in  a 
series  of  three  which  present  results 
obtained  from  the  initial  phase  of  a 
program  entitled  "Radar  Target  Char- 
acterization via  Complex  Frequency 
Resonances."  Our  objective  was  to 
assess  the  feasibility  of  classifying 
or  identifying  (C/I)  ship-like  targets 


by  use  of  radar  returns  relative  to 
the  complex  resonances  that  char- 
acterize them. 

Present  test  data  do  not  confirm 
the  feasibility  but  are  highly  encour- 
aging. The  basic  tools  required  for 
this  application,  such  as  the  pole- 
extraction  (frequency-domain  Prony) 
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technique,  and  the  predictor- 
correlator,  have  been  developed  and 
sufficiently  tested  to  show  that  the 
fundamental  concept  is  still  sound 
and  potentially  viable  for  the 
intended  application.  Significantly, 
workable  solutions  and  approaches 
appear  to  be  available  for  resolving 
questions  that  directly  concern  the 
C/I  application  (as  outlined  below) . 

Major  accomplishments  of  the  study 
include : 

• Development  and  extensive  test- 
ing against  both  analytically 
specified  data  and  computed  EM 
spectra  of  the  frequency- 
domain-Prony  technique. 

• Extension  of  the  above  for 
application  to  magnitude-only 
data. 

• Development  and  application  of 
a f requency-domain-predictor 
correlator  for  C/I. 

• Demonstration  that  iteratively 
processed  data  permits  more 
accurate  determination  of  real 
poles  as  well  as  separation  of 
curve  fitting  and  noise  poles 
(Appendix  A) . 


Significant  results,  observations 
and  conclusions  include: 

• Up  to  20  pole  pairs  can  be 
extracted  from  spectral  data. 

• Valid  poles  can  be  obtained 
from  band-limited-type  data. 

• Frequency-domain  Prony  is  use- 
ful for  interpolating  discretely 
sampled  data. 

• A fairly  wide,  dynamic  range  in 

3 

residue  amplitude  (>10  :1)  and 
damping  constant  (_>10:1)  can  be 
handled . 

• Near-  and  far-radiated  fields, 
scattered  fields,  and  surface 
current  and  charge  can  be  used 
to  extract  the  poles. 

• Comparison  of  measured  and  pre- 
dicted results  indicates  the 
most  useful  frequencies  for 
discriminating  between  targets. 

• When  a predictor-correlator 
approach  is  used,  direct  pole 
extraction  is  not  necessary. 
However,  it  might  be  preferable 
when  only  a few  poles  are  being 
excited  for  any  angle  of 
incidence . 


Recommendations 


We  recommend  the  following  areas 
for  further  study: 

• Continue  evaluation  of  frequency- 
domain-Prony  pole  extraction  as  it 


is  affected  by  noise,  source  var- 
iations, band  limited  data,  etc. 

• Continue  evaluation  of  predictor- 
correlator  strategies. 
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• Study  the  information  content 
of  scattered  data  from  the 
viewpoint  of  optimizing  dis- 
criminability  among  targets  of 
a given  set. 

• Examine  the  utility  of  pole-set 
data  and  the  predictor- 
correlator  for  developing  a 
similarity  index  between 
targets. 

• Test  the  pole  extraction  pro- 
cedure against  experimental 
data  by  using  the  iterative 
moving-window  scheme. 

• Investigate  the  possibility  for 
determining  the  information 
rank  of  spectral  data  and  how 
that  might  be  used  to  make  C/I 
more  effective. 

• Compare  the  use  of  poles  in  a 
predictor-correlator  sense  as 
was  done  here,  or  alternatively 
by  extracting  poles  from  mea- 
sured data  for  direct  compari- 
son with  a library  pole  set 
via,  for  example,  an  optical 
correlator. 

• Quantify  the  information  con- 
tent of  frequency,  transient 

* 

and  angular  data. 

* 

Information  content  might  be  quan- 
titatively defined  as  the  data-rank 
times  the  number  of  parameters-per- 
rank  times  the  number  of  bits-per- 
parameters,  for  a specified  pole  set. 
For  sampled  data,  we  might  use  the 
number  of  nonredundant  input  samples 
times  the  number  of  bits  known-per- 
sample . 


• Examine  the  feasibility  for 
direct  inversion  of  pattern 
data  by  using  a modified-Prony 
technique  to  determine  radiating 
source  locations  in  space. 

In  closing,  it  is  worthwhile  to 
comment  generally  on  target  C/I  and 
particularly  on  pole-descriptions  in 
terms  of  their  relationship  to  infor- 
mation theory.  First,  let  us  con- 
sider information  in  a more  general  way 
than  it  is  in  communications  theory, 
as  "knowledge  acquired  from  data." 

How  the  data  is  processed  to  develop 
that  information  and  how  the  informa- 
tion is  subsequently  used  provide  the 
crux  of  the  C/I  problem. 

We  have  seen  that  EM  spectrum  or 
waveform  can  be  analytically  expressed 
in  simple  terms  of  poles  and  residues. 
It  may  be  conjectured,  indeed,  that 
the  pole-residue  description  repre- 
sents the  irreducible,  minimum- 
parameter  description  of  the  spectrum 
or  waveform,  and  hence  the  EM  char- 
acteristics of  the  object  to  which 
they  belong.  In  this  sense,  it  is 
interesting  to  speculate  that  a given 
spectrum  or  waveform  possesses  a rank 
or  dimensionality  associated  with  the 
number  of  parameters  involved  in  the 
pole-residue  description.  If  this  is 
true,  and  the  rank  or  the  equivalent 

information  content  of  the  data  under 

8 9 

consideration  can  be  established,  ’ 
the  implications  for  the  C/I  problem 
could  be  profound.  For  example,  the 
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EM  rank  of  the  data  would  determine 
how  many  poles  and  residues  it  con- 
tains, or  how  many  features  it  could 
yield  in  a more  classical  approach. 

In  either  case,  its  EM  rank  imposes 
an  upper  limit  on  the  information 
that  can  be  extracted  from  the  data. 

As  we  have  seen,  noise  can  reduce  the 
rank  of  the  data,  as  demonstrated  by 
the  decrease  in  number  of  poles  that 
can  be  extracted  from  it. 

However,  if  the  problem  is  one  of 
imaging  rather  than  merely  C/I,  the 
EM  rank  must  then  equal  or  exceed  the 
object's  geometrical  rank  for  the 
goal  to  achieve  an  image  that  is  con- 
gruent with  the  actual  geometry. 

Without  restrictions  on  the  geometry, 
the  required  information  may  be  una- 
vailable from  finite  data,  i.e.,  the 
inverse  problem  could  not  be  uniquely 
solved,^'  but  the  EM  and  geometrical 
ranks  should  not  generally  be  expected 
to  be  equal.  In  any  case,  we  intu- 
itively expect  that  if  C/I  is  the 
goal  as  opposed  to  imaging,  the 
information  needs  should  be  less 
demanding.  In  other  words,  if  we 
have  some  idea  of  what  we're  looking 
for,  the  set  of  objects  of  concern 
and  the  required  information  are 
reduced . 

Thus,  the  concept  of  information 
content  in  data  would  seem  to  be  a 
central  one  in  C/l.  The  pole-residue 
description  provides  a characterization 
that  is  interpretable  in  an 
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information-content  sense,  analytically 
convenient,  which  is  why  it  is  espe- 
cially useful,  besides  being  physi- 
cally meaningful.  Because  the  infor- 
mation of  the  time-and  frequency- 
domains  is  exactly  equivalent,  they 
should  be  viewed  as  complimentary 
expressions.  The  coice  of  one  or  the 
other  should  depend  on  the  best  method 
for  acquiring  and  processing  the  data. 

Because  most  data  available  to  use 
in  the  real  world  are  obtained  in  the 
frequency  domain,  transient  waveforms 
can  be  usually  obtained  only  by  trans- 

•k 

forming  frequency-domain  data.  This 
fact  alone  favors  the  frequency-domain 
technique  for  subsequent  processing. 

In  addition,  the  frequency-domain 
method  for  prediction  and  correlation 
is  straightforward  and  makes  it  pos- 
sible to  quantitatively  determine  fre- 
quency intervals  most  useful  for  C/I. 

In  summary,  we  have  1)  developed 
and  validated  a technique  for 
frequency-domain  pole  extraction; 

2)  demonstrated  an  iterative  pro- 
cedure for  reducing  both  the  pole-set 
ambiguity  due  to  noise  and  parameter 
variations;  and  3)  shown  the  potential 
of  a frequency-domain  predictor- 
correlator  for  C/I.  Much  work  remains, 
but  progress  to  date  provides 
justification  to  continue. 

k 

At  best,  the  transform  to  the  time 
domain  can  only  preserve  the  infor- 
mation content,  and  may  actually 
decrease  it. 
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Ship-model  1 


180  m 


Ship-model  2 


30  m 


-90  m- 


Ship-model  3 


Fig.  1.  Stylized  ship  models  used  to  test  the  identification 
routines . 
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ig.  2a.  Backscattering  spectra  of  ship-model  1 for  30°  and  60°  incidence. 


30°  incidence 


60°  incidence 


Fig.  2b. 
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Backscattering  spectra  of  ship-model  2 
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for  30°  and  60° 


incidence . 
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Fig.  3a.  This  series  of  figures  shows  the  poles  and  residues  set  found  when 
the  frequency-domain  Prony  is  applied  to  the  backscattering  data  for  the 
ship-models  shown  in  Fig.  1 and  spectra  shown  in  Fig.  2.  The  data  are 
presented  in  three  separate  formats:  1)  a Cartesian  plot  of  all  left-hand 
plane  poles;  2)  an  isometric  plot  of  the  "good"  poles  in  the  second 
quadrant;  and  3)  tabulation  of  the  poles  and  residues.  The  data  shown 
here  are  for  ship-model  1 at  30°  incidence. 
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0 -7  606360E  *05 

9 -6  72I090£*O5 

13  -6  ’? I 893E  *05 

11  -6  5309606*05 

1?  -0  5S896CE  *35 

13  - 66-  3656  • 05 

- •'  66—  I 06  *05 

15  -0  60-3-6*05 

1 6 8 60— 7*6  *05 

1 7 2-60876*06 

18  -I  2-68906  *06 

!9  32907  J£. 06 

20  ; 3290756*06 

2 1 8 096-0  JE. 05 

2?  0 0965 1 26  • 05 

23  I 210077*. 06 

2-  -1  21 00706  *06 

25  0 60733-6*05 

26  -0  607J36E  *05 

27  ! 799!86€»06 

29  079',6**E»06 

30  I 07976V  *06 

31  -l  5009706  *C6 

12  . 5009706  *06 


7 5",6C88E  *05 
' 5"'6088E  ♦ 05 
: -236506  *06 
-2365C6  »06 
9-68866 ♦ 35 
- 0-68866 • 36 

2 "’-6950E  *06 
-2  ■’-6953E  *06 

3 337697E.36 
-3  30759',E*06 

- I 07602E  *06 
--  '.0-602C06 

- 656  - 96  *06 
--  656-  8E  *06 

- ■’37-0OE  *06 

- ’37--9E»06 
5 5-  3 1 86C  • 06 

•5  5-31056  *06 

5 07 1 92’E  *06 
-5  0-,:927£.O6 

6 6392726*06 

6 6392726  *06 

7 -5920 IE  *06 
-7  -592826 *06 

0 2 ! 570-6  *06 
0 21570-6  *06 

0 776052E  *06 
■«  ' ’6052E  *06 
9 207-096*06 
9 207-096  *06 

1 0II38IE»0’ 
-I  01  I 30IE  »07 


Pole  data 


Residue  data 


1 - 02900-6*01 

2 --  . 029396E  *01 

3 2 60707-6*03 

- 2 607057E *03 

5 3 62 1 9 1 3C  • 02 

6 3 62 ! 695E • 02 

7 I . 16-8806  *0- 

0 I 16-0796*0- 

9 -2 .-6762-6  *0- 

10  2 -676296  *0- 

11  - 1 6822 1 E *03 

12  - 1 60  325E  *03 

13  6 502027E  *02 

1-  6 5 1 5082E • 02 

15  -7  555-5 IE *03 

16  -7  5569 3-E  *03 

17  i 3O3-07E  *0- 

18  1.30  3 - 2 36  * 0- 

19  - 9085766*0- 

20  - 908525C  *0- 

21  -5  |7?008E*0- 

22  -5  I 720326  * 0- 

23  9 8005926*0- 

2-  9 0O0572E  *0- 

25  6 8636- 16*0- 

26  6 06S635E  *0- 

27  3 073I67E  *0** 

28  3 07?937t»0- 

29  8 74- |26E  *0- 

30  8 7—  1 53E  ♦ 0- 

31  -6  -7 1 335C  *0- 

32  -6  -7I-IV0- 


- I 5-  3 : 35E  • C 3 

1 5-31 28E  *03 
l . 02 16C8C  * G- 

-I . 021 609E • 0- 
7 330239E • C 1 
-7 . 330086E • 0 : 
-3  929085E  *03 
3 9298 3^1  • C 3 
3 60285-6*0- 
-3  6820-2E • 0- 
-2. 3518696*0- 

2 35 1 08OE • . * 
7 099-056*02 

-7  89-0-8E  • Pc 
-I  1-02326*0- 
I . l-0l-7£  *o- 
I 2790701'  *06 

- I 279083£  *05 
-I  670 I20E  *05 

1 670130f*05 

2 1 9070  IE  *0- 

2 1 9069 IE  *0- 

3 0256—1*0- 
3 0256206*0- 
- 6830— 6*0- 

-- .683I66T *0- 
-3  27  I 7881 *0- 
3 27183-1*0- 
-2  298097{  *0- 
2 29820-E  *0- 
9 38082  33  *0- 
-9  38965-f • 0- 


rijByf  *\  f T ;HC  values 


CUBVt  nrtlNO  VALUES 


33  ' ’-O07JE.O7 

3-  l 36*37 

35  5608  ’-f  •07 

36  56087-6*07 

17  O07g^gf  *Q7 

3H  067620C *07 

39  I - I5920E  *07 

-i  *159206*07 


3 93-0216*06 
- 3 93-0216  *36 

■ 

SOT*  £*M 
- 

6 2333I7£.05 
9 28  3 3306  • 06 
9 2833306*06 


33  l 8-66-26*05 

3-  l 0-665-f  *05 

35  8 0211—  E *05 

36  8 82 1 059E  *05 

37  l . I - 325 1 E *05 

30  I 1-327  |E  *05 

39  -0  15605-E  *05 

-0  H 15602-6*05 


-2  583236  *05 
2 1 50  302E  • 05 
l 3-8089C  *05 
1 3-8C96E  *05 
. 905225C • 0- 
-I  90-9506*0- 
0 5-6-966 • 05 
0 5-6-656*05 


Fig.  3b.  Same  as  Fig.  3a,  except  for  ship-model  1 at  60°  incidence. 
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X 


N 


h 
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REAL  part 

Calculated  pole  positions 


P*B* 


;-ag 


»E*w 


PART 


l mao  par* 


-2  66-5621 *;6 
-2  -T 

-1  626  36  If • ' 5 
-I  6263-31*05 
-5  '*‘*792 IE *06  2 

-5.— 791 9!  *05  -2 

--  —569-1  *06  2 

--  '•‘•668: E *06  -2 

-9  -5268-1*05  3 

-9  -5268-1  *05  -3 

-6  1657551 *05  - 

-6  1657561 *05  -- 

-0  66199-1*05  - 

-9  66 1 99 E • 06  ->* 

85-9901*06  6 

i 93-9911 *06  -5 

6993'  Sf  • 36  6 

-7  69936**E  *06  -6 

-I  29-2691*06  6 

-»  29-2681 *06  -6 

•7. 7 1 09991  *06  7, 

-7 . 7 i 0986E  *06  -7 

*1.1 6062'*E  • 06  8 

- I 160620C  *06  -8 

■.  \ 8 

: 3-?79ll  *06  -8 

-I  7Q678 If  *06  9 

9 

3 9101-81 *05  J 

30  -3  8101 60E » 06  -I 


88' 39-1 *05 
88' 386E • 35 
2'*699"'E  *06 
2>*699'E  *06 
10-3561  *36 
8-  3661  ♦ 06 
-603 . 8(  • 36 
-603.01*36 
-807921 *06 
'•807921  *06 
l C3607£  *06 
l 0368'E  *06 
7-26591  *06 
’-26591  *06 
55 36: -1*06 
6636:  ■•£•06 
763796£*06 
763"’96€  *36 
76l636f  *06 

Pole  data 

-216201*06 
1 0262"’E  *06 
1 02627E  *06 
8968'*8£*06 
0968**  ’E  • 06 
006-01  *06 
l 086**  8f  >06 
0006601*07 
0006601*07 


Residue  data 


: 03611 91 *03 

2 1.8360991*03 

3 6 96367 IE *02 

- 6 96371 OE *02 

6 -1 . 166735€*09 

6 -l  1567331*0- 

7 1 .50—011*03 

0 l . 60'*'*  l 3E  *03 

9 -.367*.  991  *02 

10  -.36 76 391 *02 

11  2 6070091*03 

12  2 6070161 *03 

13  - 1 230681 ‘O'* 

1-  --  1 230661 *0- 

16  2 260682E  *0- 

16  2 2606921*0- 

17  3 -966 1 7E  *03 

18  3 -9661 3E*03 

19  3 66— 69E*03 

20  3 66-626C  *03 

21  3 1-O076E  *03 

22  3.1-09301*03 

23  1 93-6671*03 

2-  I 93— 851  *03 

26  -2  0661661*0- 

26  -2  05601 31*0- 

27  --.0720931*0- 

28  - 8 9861 *0- 

29  I -910881*03 

30  I -910—1*03 


-7  1661 30E»03 

7 1560291 *03 

5.-09-201*03 

-5  -09-361*03 

6 "’7  32911  *03 
-6  7732791*03 
-6  - 7836-E  *03 

6 -703671 *03 

-7  709967E.O2 

7.7000—1*02 

-l  I 159261 *0- 

1 : 15927E  *0- 

2 326 1291*0- 
-2  326 1261*0- 
-l  2133131*0- 

1 2133231*0- 
-I  1731 371*0- 

1 . 1 73 1 32E  *0- 

- 9069l5£*03 
--  9069121*03 
-I  26-3101.0- 

l 26-3031 *0- 

- 05',689t  *03 

- 0579971*03 
-1 .96-7031 *03 

1 966-681 *03 
-3.26991 11*0- 

3 26966 11*0- 

-2. 2703291*03 

2.277'331*03 


C <»v€  f I • ' NC.  » At  UCS 


CUfiVE  FITTING  VALUES 


31  9 5790961 *06 

12  9 5790961 *06 

3»  - 68—0-1*00 

3-  ••  68— 0-C  *06 

36  1 9211  61 *06 

R • ’ ■ • • 

37  ! -76-911*0" 

38  1 -76-9lt*C7 

JR  1 0122671*06 

-0  3 0122671*06 


5 291 1701*06 
-5  291  I 70E  *06 
3 1661  131*06 
3 160  31*06 

: 2977091 *06 
1 29"C9f  • 6 
' 280-261 *06 
...  . 

■ m ■ • 

-9  380 : 701*06 


31  . 665-1 11*0- 

32  -1  6566651 *0- 

33  -2  60386-1*03 

3-  2 6039801*03 

36  2 5063371*03 

36  2 5062-61 *03 

37  1 3660761 *05 

30  1 3669101*05 

39  5 3100901*03 

-0  5 3102671*03 


. L.t-J  1 

--  3— 6061*0- 
9 6-093-3  *03 
-5  6-C'38E*03 
-2  -0-0161*03 
2 . -0500-t  *03 
9 —39701*02 
9 — 37221*02 
-3  5353901*03 
3 53,939C  *03 


Fig.  3c.  Same  as  Fig.  3a,  except  for  ship-model  2 at  30°  incidence. 
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l . OC  ► • 

u . 75  * 

0 50*-  , * 

o 06*  . • 

0 oo  ► . 

-o  95*-  '•  . J 

-0  60  * 

-0.75f- 

. 

1 05 1 *-  - * — - . . l J 
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real  p art 

Calculated  pole  positions 


Isometric  plot  of  calculated  poles 


***G 


p*o* 


»C*L  PART 


■£  89:  0 ; 9C *C5 

9 £ 09iOI0E*O6 

3 . ■’?3-S8f*^5 

8 . '63-  591  * C5 

6 6 66  996f • 35 

b -6  55l909f*03 

-6  0 33',39€  *05 
0 -6  C33’39€*S5 

9 -9  3??9;-E*05 

tO  9 0999-1*06 

II  -6  986771*05 

:9  6 <9M77c*05 

J 980  I 39E  *06 

|8  98C:?9£*36 

i 6 • 968-E*06 

.6  '.86881*06 

. • 6 9906 1 9C *05 

10  -6  99060 IC  *06 

;9  : 33^6636*06 

90  -i  S38668i*06 

91  -1  769977E*06 

99  1 7W9fft  *06 

93  1 9369; 3E  *06 

08  : 93691  If *06 

96  ! 96  393*»t  *06 

96  96  39 J6€  *06 

9’  7 9696 3'*f  *06 

98  ' 9696 36f  *06 

99  ' '*6*t3“t»06 

10  I -6-!3-»f*06 


7 88-8'0£*O6 

' 88-e'8E*36 

. 963063C  *06 
963063E  *06 
9 :0'**'6E  *06 
-9  I8^«76£*06 
9 ‘*-9996£*C6 
9 8-009V  • 06 
J - 89801  T6 
3 - 89801  • 06 
M - - 
- :0C003C*06 

8 7h  39**h£  *06 
-•*  78  S088E  *06 

5 5899**h£  *06 
* 68998-1  *06 

6 75907w£ .36 

6 769C7««E.36 



-6  7800001*06 

7 9 1 9986f  • 06 
-7  - 9986C  *06 

8 I90697E  *06 
8 1906971  *06 

8 9806081  • 06 

B 98969xE*06 

9 07-3981  *06 

I 0 1 t)7"'8t  *07 
OlO-’70£  .07 


Pole  data 


Residue  data 


I 6 707|3**E*09 

9 5 . 707 1731 » 00 

3 -0083181*03 

>*  - 99u  39 1 E ♦ 0 3 

5 -9  85J088E *03 

6 -9  -55036C  *03 

7 6 C78999E.03 

0 6 0 "’66 ! 0E  ♦ 0 3 

9 -l  7867801*08 

10  •!  7467V.0E.O1* 

11  I 9603ClE*0‘* 

19  ; 960309E  *0** 

13  -6  1 6 1 599E ♦ 09 

l"  6 16I697E.09 

16  -I  099308€  *05 

■ 6 . 09930BE  *06 

17  : 93083V  *38 

IS  I 9308511*08 

19  -I  0088071  *05 

9C  -1  00800V  *06 

91  9 09 1 769E  • 0** 

££  9 . 09 1 76 1 E • 0** 

93  * 3698981  *03 

0“  8 3697991  *03 

96  6 7501311*08 

96  -6  75799OC  *0** 

97  1 670CI6E  *0** 

08  1 669966C  • O'* 

99  1 08- 3391 *05 

30  1 08-3301*05 


-1  970-061 *C! 
1 978-0-1  *03 
3 988-*76c*C3 
-3.088’'9E*C3 
-8  -555331  *03 
8 -555661 *03 
-1  969360C • 0* 
1 969 368E  ♦ 0- 
- 0693891 *03 
-8  969376E  *03 
-8 .39330  IE *08 
8 3933031  *08 

0 O07098E  *0 1 
-9  08091 !E*01 

8 6898"' it  *08 
-8  68083-1*0- 
-5  679806C«08 
6 670386? -08 

1 0080971*0- 

62691  • 8 
6 1 30 1 69C  • 08 

6 I30I68£.08 

-0  808669f • 

0 808687C. 

1 7968 VHE* 

- I 7968  1 St  . 

-9. 39996 IE* 

0 399666C • 08 
3 09301 7f  *08 
3 09308  IE  *C*8 


VALUES 


'U»vt  f MTING  VALUl  6 


3i  6 69  388 7£  *06 

39  6 6- 388  T*  06 

JJ  6 668P08f  *06 

38  6 568008f  *06 

36  3 800  3- Cl  *06 

36  6 ’8i ; V36 

37  

38  6 781  1 6C  • 06 

39  1 061  T60f  *08 

80  06  1 760f  » 08 


6 966660E • 06 
6 966*>60C*  6 
0 '038601  *06 
0 7c’3060C*OS 
6 C9-598C  - 18 



60690V  39 



3 807660C*06 
3 8Q7660E  *06 


31  6 0606361*03 

30  -6  0606181*03 

33  I 7861901*03 

38  1 7^*^071*03 

36  8 706  3991  *00 

36  0 980w06€  *03 

37  8 385  396E *C3 

3B  0 0-00831*03 

39  1 387 3801* 06 

8fl  | 3878  781.06 


1 6089931  *03 
1 6093601*03 
8 3 : 7-.50C  .03 
8 317506E  *03 
•I  30068V  -09 
6 9107681.03 
6 61 18691 -00 
-6  9l8080t*03 
5.912  • I 

6 9' 96961 *06 


Fig.  3d. 


Same  as  Fig.  3a,  except  for  ship-model  2 at  60°  incidence. 
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BE 


e»p* 


*1». 


P*s- 


“*G 


! 1 ‘•3-90',E*2  3 

2 : W5^917E.03 

3 1 ; 989  3-*€  *04 

••  98894t*.> 

9 -2  S6e686E*:«. 

6 -2  S62680E*04 

7 -0  ‘♦0*0296*03 

0 -0  4040:91*03 

9 6 9022.  c+03 

10  6 9829491 *03 

11  2 '946801*34 

I 2 2 1946071*04 

: 3 -2  .837921*04 

;*•  -2  i03‘’8^E»C'* 

19  -8  V 2 821*02 

,6  8 9! '.7671*02 

i7  7 070246i*c3 

.0  7 0-03491 *03 

19  1 2.4-891*04 

20  : 2 : 4-991  • S'* 

21  2 9246401*04 

22  2 9246801 *04 

23  2 0427831*0  3 

2'*  2 0428981  *03 

28  -2  221  3i-4£*03 

26  -2  221 0911*03 

27  9 308 1 991  *03 

28  9 3082831*03 

29  -8  33l093E*04 

30  8 33!CB:£*0« 

3:  -•*  m3':,*E*09 

32  •«  i-*3-2?E*08 


-h  3299"E»03 
•*  3288  38E  *03 
, 8870: ie*C“ 
»7f7flt*0H 
2901*C0E*0- 
: 2909381 • 0“ 
0 0999921*03 

0 899-33t*03 
-I  60-9361  • : - 

6:'9’».  • •• 
-I  1 96298E  • C ** 
96  I 

* 8802  *£•:- 
-•*  888c  ’1*0" 

-2  '*>.920-E*r$ 
2 •*^92C-E*0 J 
-2  -269  391*0- 

2 ■’268  39C  • O'* 
, 672397T*> 

-I  6-23911*0- 
9 88632  JC  *C  * 
9 886  36-1  *0  3 

1 20-382€*0- 
-1  204  34-t.r-* 

4 269602E  • C ‘ 
- 2S9’88E • : 3 

3 9866  3-1* 

3 96-0991 *02 
8 9942011  O- 
-6  99-3  91  • - 
3 S98-6Cf • 0m 
3 098^681 -08 


Pole  data 


Residue  data 


02- 6221*09 

2 02-6261*09 

3 -4  789:  161*08 

4 - -89: 081  *08 

8 -2  9:94601*06 

& -2  9 1 8**83E *06 

-8  882620E  *08 
0 9 8028391 *08 

9 ;08>*86E*C6 

;0  : 08" 921 *06 

l;  -6-386601*08 
12  -6  -386601*08 

3 9 99888hE*08 

!u  -g  99888*1 • 08 

-3  *0249-1*08 
!6  3 -824966*08 

3-  18991 *06 

10  -1  3- 18991 *06 

9 : 6860981*06 

20  : 6860961 *06 

21  0032891*06 

22  : 0032621-06 

23  9 —90201*06 

2-  -9  — 9Cb0£*O9 

28  8 807,891*08 

26  9 88--  1*08 

27  -J  4918631*09 

28  3 49:9931  *09 

29  -2  809  391*06 

SO  2 80*  39<  - * 

31  f 1939791 -06 

32  -6  1639-91 *06 


- 4342251.09 
--  4342251*06 

1 2-00:91*06 
2700131*06 

1 3002931*06 
3002631*06 

2 293:631-06 
-2  293:841*06 

2 -290101*06 
-2  -290081*06 

- 0392731*06 
-4  0392-31*06 

- ‘•021021*06 
-‘•.‘•021021*06 

9 1 688-01 *06 
-6  : 698-01 *06 
9 -994831*06 
9 -994831*06 

6 697-681 *06 
-6  6974681*06 

' 2373221*06 
2373221*06 

7 998-901*06 
' 996-901*06 
9 H 93- 361  *06 
S -93- 361 *06 
9 101-8-1*06 

- • • 

9 4*:eesK  *06 
9 45ieeet*06 


. 04.’40t*0- 


CjRvt  f ; 7-  i*«o  .*l  j«s 


<jBvt  ri"iNG  w^uts 


33  -4  C'97-6€*0S 

34  4 0-8808C  *C  3 

39  9 90'  991  *02 

36  9 98-4621*02 

37  3 2792  381  • O4* 

38  3 2-9  '01 *0- 

39  -2  6029841*0- 

40  : 422-1 :E*o9 


3 '628301*03 
S '62'29t*03 
■ RMt< 

-2  0-93421*0. 
-2  909: -CE  *09 
2 909: 4 1C *09 
2 344 > 201 -C: 
2 90 i 77|E -o : 


33  2 4699331*06 

- 2 4699331*06 

39  8 693-: 21*02 

36  0 6980921*02 

3’  : 4302481 *C- 

30  I 4 302481  • C ' 

39  6 0-3-661 *06 

40  l 6919901*0' 


6 23800-1*06 
6 23808-1*06 
3 3-4 !02t *06 
. ...... 

' 660931 *06 

‘ 660931*06 

3 '239:31  -04 
--  0949441-09 


Fig.  3e.  Same  as  Fig.  3a,  except  for  ship-model  3 at  30°  incidence. 
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-40’ 


:-*o 


o*«r 


P*0’ 


MAG  PA»’ 


-8  083  O' 6*05 
■a  0831016*35 

0 9O6-50T  *P- 
-9  9063396*:)- 

... 

-5  B3i‘:et.;5 
: J0C9986O6 
- 3009086 *06 

1 06-1636  *:s 
06-.  6 36*06 
9889556*05 

• 5889656* 05 

5 35068  36*05 

6 95868-T  05 
-0  'TOC’t  05 
•0  •O05C-EO*. 

0 05-3186 *06 
• . • ... 

1 ,o?m-’:e.o6 
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Fig.  3f.  Same  as  Fig.  3a,  except  for  ship-model  3 at  60°  incidence 
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Fig.  4a.  This  series  of  figures  shows  the  results  obtained  when  the 
predictor-correlator  algorithm  is  used  to  identify  a scattered  re- 
sponse by  three  sets  of  stored  library  poles.  The  backscattered  data 
shown  in  this  figure  are  for  ship-model  1 at  60°.  The  three  sets  of 
library  poles  used  in  the  identification  algorithm  are  for  ship- 
models  1,  2,  and  3 at  60°.  This  figure  shows  the  magnitude  data  of 
all  points  used  for  both  the  calculation  of  residues  and  the  sub- 
sequent identification  calculation. 
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Fig.  4b.  The  magnitude  plot  of  the  scattered  response  shown  in  Fig.  4a 
for  the  points  where  the  residues  will  be  calculated  is  shown. 
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Fig.  Ad.  The  magnitude  plot  of  the  calculated  identification  response 
for  ship-model  j.  is  shown.  This  response  is  found  by  using  the  stored 
library  poles  for  ship-model  1 and  the  data  in  Fig.  Ab  to  calculate  a 
set  of  residues.  This  set  of  residues  and  the  stored  pole  set  are 
then  used  to  calculate  a response  function  at  the  identification  points. 
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Fig.  Ae.  The  magnitude  plot  of  the  difference  between  the  measured  and 
calculated  responses  at  the  identification  points  is  shown.  This 
function  is  found  by  subtracting  the  response  used  to  find  Fig.  4d 
from  the  response  used  to  find  Fig.  Ac.  The  identification  number 
[Eq.  (4a)]  with  ship-model  1 is  equal  to  0.0. 
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Fig.  4f.  The  magnitude  plot  of  the  calculated  identification  response 
for  ship-model  2 is  shown.  This  response  is  found  by  using  the  stored 
library  poles  for  ship-model  2 and  the  data  in  Fig.  4b  to  calculate  a 
set  of  residues.  This  set  of  residues  and  the  stored  pole  set  are 
then  used  to  calculate  a response  function  at  the  identification  points. 
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Fig.  4g.  The  magnitude  plot  of  the  difference  between  the  measured  and 
calculated  responses  at  the  identification  points  is  shown.  This 
function  is  found  by  subtracting  the  response  used  to  find  Fig.  4f 
from  the  response  used  to  find  Fig.  4c.  The  identification  number 
[Eq.  (4a)]  with  ship-model  2 is  equal  to  0.156. 
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Fig.  4h.  The  magnitude  plot  of  the  calculated  identification  response 
for  ship-model  3 is  shown.  This  response  is  found  by  using  the  stored 
library  poles  for  ship-model  3 and  the  data  in  Fig.  4b  to  calculate  a 
set  of  residues.  This  set  of  residues  and  the  stored  pole  set  are 
then  used  to  calculate  a response  function  at  the  identification  points. 
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Fig.  4i.  The  magnitude  plot  of  the  difference  between  the  measured  and 
calculated  responses  at  the  identification  points  is  shown.  This 
function  is  found  by  subtracting  the  response  used  to  find  Fig.  4h 
from  the  response  used  to  find  Fig.  4c.  The  identification  number 
[Eq.  (4a)]  with  ship-model  3 is  equal  to  0.122. 
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Appendix  A 

Repeated  Processing  Cluster  Plots 


The  presence  of  curve  fitting  and/or  noise  poles  due  to  overdetermining 
the  system  is  difficult  to  detect  in  a single  calculation  from  a given  data 
sample.  If  additional  sample  sequences  are  generated  from  the  same  waveform 
(time-domain  Prony)  or  spectrum  (frequency-domain  Prony) , it  appears  feasible 
to  differentiate  between  real  poles  associated  with  the  object  and  anamolous 
poles  generated  by  the  process. 

The  additional  sequences  might  arise  in  several  ways.  Perhaps  the  sim- 
plest in  the  frequence  domain  is  the  addition  to  random  numerical  noise  to 
the  original  date  to  generate  a new  sequence.  The  invalid  poles  presumably 
would  move  in  response  to  the  noise  more  than  would  the  real  poles,  thus  pro- 
viding a means  for  separating  them.  Repetition  of  this  process  several  times 
should  produce  clusters  of  real  poles  and  expedite  locating  them.  This  tech- 
nique could  tend  to  degrade  the  overall  process  however,  and  would  have  to 
be  exercised  with  care. 

An  alternate  procedure  more  analogous  to  what  has  worked  in  the  time 
domain  would  be  to  develop  the  additional  samples  by  moving  the  frequency 
interval  over  which  the  poles  are  calculated.  Again,  by  repeating  this  proc- 
ess some  number  of  times,  we  can  identify  the  valid  poles  falling  in  the  fre- 
quency range  overlapped  by  all  intervals. 

A third  approach  possible  in  the  frequency  domain  would  be  changing  the 
angle  of  incidence  and/or  viewing  angle.  An  advantage  of  the  former  is  that 
poles  not  appreciably  excited  at  the  first  angle  of  incidence  chosen  would  be 
less  likely  to  be  missed.  The  clustered  poles  that  result  should  be  repre- 
sentative of  a truly  aspect-insensitive  pole  set,  and  thus  minimize  the  cor- 
relator uncertainty  discussed  relative  to  Table  1. 

Examples  of  some  of  these  possibilities  are  shown  in  Figs.  A-l  to  A-3. 
In  Fig.  A-l,  we  present  the  cluster  plots  that  occur  from  a fixed  bandwidth 
window  being  systematically  shifted  through  the  spectrum  being  processed. 

The  target  in  this  case  was  ship-model  3 for  0 = 60°,  with  data  available  at 
0.125  MHz  intervals  from  0-10  MHz.  A bandwidth  of  5 MHz  was  employed  for  the 
pole  calculation,  resulting  in  40  complex  data  points  and  20  pole  pairs. 
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The  poles  were  derived  from  ten  separate  calculations  by  using  the  fre- 
quency intervals  0.125(i  +1)  to  5 + 0.125i  MHz,  i = 0,...,9.  Figure  A-l 
shows  the  pole  sets  obtained,  where  the  pole  locations  are  designated  by 
A,...,  J corresponding  to  i = 0,  i=l,....  Most  of  the  poles  actually  present 
in  the  spectrum  are  well  defined  by  clusters  whose  size  is  that  of  the  indi- 
vidual letters.  There  is  greater  variance  in  some  locations,  such  as  for 
W > 5.  Evidently,  they  are  initially  outside  the  window  of  data.  Curve- 
fitting poles  usually  seem  to  scatter  rather  randomly  about  the  s-plane. 

An  example  of  varying  the  incidence  angle,  again  for  ship-model  3,  is 
shown  in  Fig.  A-2.  A bandwidth  of  0-5  MHz  was  employed,  with  40  frequency 
samples  at  0.125  MHz  intervals  and  20  pole  pairs  being  extracted.  Incidence 
angles  of  9 = 25°  to  65°  in  5°  steps  were  employed,  and  the  corresponding  pole 
locations  designated  by  letters  A,...,J.  For  those  poles  whose  location  are 
given  by  letter-size  clusters,  the  agreement  with  Fig.  A-l  is  essentially 
exact.  Curve-fitting  poles  are  made  rather  obvious  by  the  randomness  of  their 
locations . 

As  the  last  example  of  cluster  plotting  given  here,  we  present  in 
Fig.  A-3,  the  results  of  using  a shif ted-f requency  window  for  data  contamina- 
ted by  noise  (ship-model  3 at  0 = 60°).  Cluster  plots  are  shown  for  three 
cases  of  additive,  white  noise  having  a maximum  value  normalized  to  the  peak 
of  the  spectal  response  of  ~1.25%,  3.75%  and  11.25%,  respectively.  With 
increasing  noise  levels,  the  locations  of  valid  poles  become  more  uncertain, 
but  the  effect  depends  strongly  upon  the  residue  magnitude  associated  with  a 
given  pole.  Naturally  enough,  those  poles  with  the  largest  residues  are  least 
sensitive  to  a given  amount  of  noise,  while  poles  having  residue  values  near 
the  noise  level  may  disappear.  In  other  words,  when  the  noise  is  comparable 
to  the  magnitude  at  which  a given  resonance  is  excited,  that  resonance  becomes 
unrecoverable  from  the  data.  The  rate  at  which  information  is  lost  from  a 
spectrum  or  waveform  as  noise  increases  would  be  useful  to  establish  in  a more 
quantitative  sense.  For  the  shif ted-window  approach  to  be  successful,  uncor- 
related noise  or  a valid  noise  model  is  required.  This  is  also  true  of  other 
methods  for  processing  noisy  data.  The  repetitive  processing  of  such  data 
effectively  increases  the  amount  of  information  available,  and  thus  improves 
the  numerical  results  that  can  be  obtained. 
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Fig.  A-la.  A frequency  cluster  plot  for  ship-model  3 at  0 = 60°. 
This  plot  is  produced  by  using  a 5-MHz  frequency  window,  10 
different  times  on  the  spectrum  data  for  0-10  MHz.  The  poles 
of  each  consecutive  run  are  labeled  A to  J. 
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Fig.  A- 2a.  An  angular  cluster  plot  for  ship-model  3.  This  plot 
is  produced  using  the  frequency  range  of  0-5  MHz  to  obtain  the 
20  pole  pairs  but  for  incidence  angles  from  25°  to  65°.  Ten 
different  runs  were  thus  performed  and  they  are  lettered  from  A 
to  J on  the  plot. 
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Fig.  A-2b.  Same  as  Fig.  A-2a,  except  this  figure  shows  only  the 
second  quadrant  on  an  enlarged  scale. 
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Fig.  A- 3a.  A frequency  cluster  plot  for  ship-model  3 at  60°  when 
the  data  are  contaminated  hy  white  noise  having  a maximum  value 
normalized  to  1.25%  of  the  spectrum's  peak. 
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Fig.  A-3b.  Same  as  Fig.  A-3a,  except  this  figure  shows  only  the 
second  quadrant. 
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Fig.  A-3c.  A frequency-cluster  plot  for  ship-model  3 at  60°  when 
the  data  are  contaminated  by  white  noise  having  a maximum  value 
normalized  to  3.75%  of  the  spectrum's  peak. 
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Fig.  A-3e.  A frequency  cluster  plot  for  ship-model  3 at  60°  when 
the  data  are  contaminated  by  white  noise  having  a maximum  value 
normalized  to  11.25%  of  the  spectrum's  peak. 
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Fig.  A-3f.  Same  as  Fig.  A-3e,  except  this  figure  shows  only  the 
second  quadrant. 
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Appendix  B 

Frequency-Domain  Predictors 

(E.  K.  Miller  and  R.  M.  Bevensee) 


A time-domain  linear  predictor,  as  expressed  by  Eq . (3),  requires  as 

data  N-constant  coefficients  a„,  , which  are  determined  by  the  system  poles 

N or 

and  N-past  samples  of  the  time-varying  waveform.  Once  the  necessary  informa- 
tion is  available  to  put  the  predictor  into  operation,  all  future  values  of 
f(t)  can  be  found  from  repeated  application  of  the  predictor.  As  used  in 
Ref.  6 for  object  C/1  however,  only  measured  data  are  inserted  into  the  pre- 
dictor while  the  computed  value  at  the  next  time  step  is  used  only  for  corre- 
lation. In  this  way  the  data  used  for  the  prediction  are  continuously  updated, 
and  computation  of  future  values  of  f(t)  involves  extrapolation  over  only  one 
time  step.  This  tends  to  give  equal  weight  to  all  the  measured  data. 

A frequency-domain  predictor's  concept  and  utility  are  not  so  obvious. 

It  is  possible  to  obtain  (at  least  implicity) , a frequency-domain  predictor 
that  operates  in  the  same  way  as  Eq . (3).  The  important  difference,  however, 
is  that  its  coefficients  depend  upon  frequency;  thus,  it  may  be  termed  a 
vary-linear  predictor  (to  borrow  from  system's  nomenclature).  Its  development 
is  indicated  below. 

First  note  that  the  time-and  frequency-domain  explicit  pole  representa- 
tions are 


N 

no  - £ 

a=l 


Srvt  & 

R e a—  — F (ju)) 
<S- 1 


2 

a=l 


R/(jto  - s ) 

a J a 


(Bl) 


In  a similar  fashion,  their  corresponding  predictor  or  implicit  repre- 
sentations take  the  form 

N N 

F(t)  = ^ ^ a(sB)f(t  - aAt)  ^=TF(jcj)  = ^ \ a(w)F(jwa) , (B2) 

a=l  a=l 

where  the  right  hand  side  of  Eq . (B2)  has  yet  to  be  derived. 

There  may  be  several  ways  to  accomplish  this.  We  choose  to  employ  the 
right-hand  side  of  Eq.  (Bl) , and  to  obtain  the  R ^ from  samples  of  F(ju^)  = M^, 
k = 1 , . . . ,N . Thus ; 
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N 

R.  = y Y M 
k ka  a 

a=l 


(B3a) 


where 


' - *_1  \cf  <J“k  - V1' 


(B3b) 


Formally,  this  leads  to 

N 


F(i“>  - I (I  VaV(J“  - 

k=l  ' a=l  ' 


(B4) 


which  becomes,  upon  interchanging  the  order  of  summation. 


N 

v- 


r N 


F(jw)  = 

a=l 


I 

k=l 


Y /(jw  - s ) 
k°</  k 


H 


(B5) 


Thus,  from  Eqs.  (B2)  and  (B5) , 

N 


Y«<“>  ' I V«“  - V 


(B6) 


k=l 


Clearly,  is  frequency  dependent,  due  to  the  ju  term  in  its  denominator. 

In  addition,  we  observe  that  since  X depends  upon  the  sequence  to,  , k = 1 N; 

k 

any  change  in  the  data  vector  M (which  thus  changes  the  to,  ) would  result  in 
= k _ k 

a new  matrix  X and  thus  a new  inverse  Y.  From  a computational  viewpoint,  the 

latter  could  be  the  most  costly  effect  of  updating  the  data  since  it  varies 

3 2 

as  N ; whereas,  a change  in  to  results  in  a N-  or  N -updating  process,  depending 

upon  whether  an  explicit  [Eqs.  (Bl)  or  (B4)  ] or  implicit  [Eqs.  (B2)  or  (B5)] 

form  is  used.  It  thus  appears  that  an  M-point  computation  using  Eq . (B2) 

3 

could  involve  from  MN  to  MN  operations,  depending  upon  whether  M is  updated 

k 

for  each  computation  or  not. 

However,  if  the  updating  is  done  in  a systematic,  stepwise  fashion,  the 
order  of  the  computation  can  be  decreased.  For  example,  suppose  the  initial 

data  vector  is  (k  = 1,...,N),  the  second  is  (k  = 2, . . . ,N  + 1)  , ..., 

(j)  k k 

(k  = j , . . . ,N  + j - 1) , where  the  superscript  denotes  the  step  number. 
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Then,  the  matrix  equivalent  of  Eq . (Bl)  at  step  j is 

,(j) 


y Y(j)  R(j)  _ M ^ 

£ xkct  Ra  - \ 


(B7) 


where  depends  implicitly  upon  j through  oo  . But  when  j -*•  j + 1 , we  note 

ka  k 

that  X is  simply  changed  by  deleting  the  top  row  (k  = j ) and  adding  a new 

old 

bottom  row  (k  = i + N - 1) . 

new 

Suppose  then  the  original  matrix  is  conceptually  depicted  as 


;U) 


x(1) 

x(1) 

...  x(1) 

11 

12 

IN 

...  x(1) 

21 

22 
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\ 

2N 

1 

1 

| 

\ 

\ 

vu> 

yd) 

V v (1) 

hi 

XN2 

• ••  NN 

and  is  solved  by  Gaussian  elimination  to  obtain  the  form 

O' 


X^' 


(1)' 

11 

(1) ' 

X(1)' 

21 

1 

X22  — 

1 

(1) 

„d) 

4XN1 


^2 


(B8) 


(B9) 


from  which  the  ^ element  can  be  found  by  back  substitution.  If  we  now 

introduce  the  k = N + 1 equation  while  dropping  the  k = 1 equation,  we  obtain 


;(2) 


vd)’ 

(1)' 

X21 

X22 

(1)’ 

(l)’ 

X31 

X32 

(1)’ 

(1)’ 

hi 

XN2 

(2) 

v(2) 

Vl.l 

XN+1 , 2 

c(lV 

33 


o 


(BIO) 


Upon  completing  the  diagonalization , we  find 


-47- 


= (2)  ' 

X 


(2)’ 

11 

(2)' 

(21 


X. 


(2) 

N+l , 1 


,(2)  ’ 
^22 


\ 


(2) 

*N+1.2 


o 


(2) 

‘N+1,N 


X 


(Bll) 


a process  which  involves  ~N  /2  operations. 

By  induction,  we  thus  conclude  that  N solutions  of  a 2N  x N system  in 

3 

terms  of  sequential  N x N subsets  will  involve  ~4N  /3  operations  in  total 
4 

rather  than  ~N  /3  which  would  result  from  not  exploiting  their  common 


equations . 

We  can  summarize  the  order  of  the  predictor  process  as  follows: 


Time  domain 
Freq.  domain 


Explicit 

Updated Not  updated 

4N3/3  N3/3 

4N3/3  N3/3 


Updated 

3 

N / 3 
4N3/3 


Implicit 


Not  updated 
3 

N / 3 
3 

N 73 
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Appendix  C 

Pole  and  Residue  Extraction  from  Magnitude  Data 


A method  to  extract  the  poles  and  residue  from  only  the  magnitude  data 
of  some  system  response  is  presented  in  this  appendix.  This  method  is  based 
on  the  assumption  that  the  response  function  is  causal  and  positive-real. 

The  positive-real  function  has  all  poles  and  zeros  in  the  left-half  plane. 

The  assumption  of  a positive-real  response  is  only  applicable  when  the  system 
response  are  from  driving-point  impedance  and  admittances.  The  numerical  pro- 
cedure will  first  be  derived  followed  by  some  results  of  computation. 

The  response  function  for  an  object  having  only  poles  and  residues  can 
be  written  as 


N 

F(s)  - £ 

a=l 


(s  - 


s ) 

a 


(Cl) 


For  the  case  of  casual  response,  the  residues,  R , and  poles,  s , occur  in 

a a 

complex-conjugate  pairs.  It  has  been  shown  in  Eqs.  (3),  (4)  and  (5)  of  a 
2 

previous  report  that  (Cl)  can  be  written  as 


F(s) 


N-l 


a=0 


(C2) 


where  A and  B are  real  numbers,  and  B = 1 . An  alternative  form  of  the  same 

a a N 

function  is 


F(s) 


N-l 

n 

a=i 


(s  - 


za^ 


N 

(s  - s ) 
a 

a=l 


(C3) 
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where  z are  the  zeros  of  the  response  and 
OL 


di  an-i 


In  Eq.  (3)  of  the  previous  report  , it  has  been  shown  that 


vr  1 


Therefore,  for  casual  systems,  the  A^.  ^ are  negative  real  numbers. 

The  magnitude  square  of  the  function  can  be  found  by  multiplying  (C3) 
by  its  complex  conjugate  to  obtain 


|F(Ju) 


, 2 a=0 


aa(iw) 


> 

a=0 

where  b = 1.  The  (joo)  has  been  substituted  for  s in  the  above  equation  since 
N 

the  frequency  response  is  known  only  along  the  jo)  axis.  The  a and  b are 

a a 

real  coefficients  since  the  A and  B were  real, 

a a 

Note  that 


d = A = a 
1 N-l  N-l 


A different  representation  for  the  magnitude  square  can  be  found  by  multi- 
plying (C3)  by  its  complex  conjugate  and  substituting  ju)  for  s.  This  new 
equation  is 


! F(  j£j) 


~di  n [ (jw  - za) + za) ] 


n[  (jca  -s  ) (ju>  + s ) ] 

a a 
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It  can  be  shown  from  the  derivation  of  (C5)  and  (C7)  that 


N-l  * 

V a (jco)2a  = -d"  FI  [(ja)  - z ) (ju  + z ); 

/h  oi  1 4 * a a 

„ a=l 

a=0 


(C8) 


and 


I 

a=0 


b^(jw)2a  = n [ (jai  - s ) ( jto  + s*) 

a ii  a a 


(C9) 


a=l 


It  is  obvious  from  the  two  above  equations  that  the  roots  of  the  two  polynomial 
equations  will  yield  the  poles  and  zeros  of  the  original  function.  To  obtain 
these  roots  requires  the  knowledge  of  the  coefficients  in  the  polynomials. 

The  coefficients  of  (C5)  can  be  found  by  writing  this  equation  in  the 

form 


N-l 


N-l 


| F ( jtd)  | 2 ]>  b^(ja))2a  + ^ aa(jio)2w  = | F ( jto)  | 2 (joj)2N . (CIO) 


a=l 


a=l 


By  evaluating  this  equation  at  2N  points  along  the  j CO  axis  the  coefficients 
can  be  expressed  as  the  following  set  of  a linear  equations. 


[A]  [C]  = [G], 


(Cll) 


where 


r -|F(jCD  )|2(j(0  )2(P_1) 


A = i 
qp 


, for  q = 1,2 2N, 

P = 1,2 N 


. 2(p-l)-2N 

(jto^)  , for  q = 1,2.  . . , 2N , 

p = N+l , . . . , 2N ; 


(C12) 
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r 


, for  q = 1 N 


and 


, for  q = N+l , . . . , 2N; 


(C13) 


F(ju>q)|2;  q = 1,2 2N . (C14) 

The  | F ( j oo  )|  are  the  values  of  the  magnitude  square  at  the  frequence  test 

point  jio  which  is  assumed  to  be  given.  The  coefficients  are  found  from 

inverting  this  set  of  linear  equation. 

Once  the  coefficients  b are  known,  the  equation 

a 


N 

S ba(jw)0t  = 0 (C15) 

a=0 

can  be  solved  for  the  unknown  roots  which  gives  the  poles.  For  causal 
responses  the  poles  of  the  response  are  in  the  left-half  plane.  It  is 
obvious  from  (C9)  that  the  roots  of  (C15)  will  be  the  original,  poles  plus  a 
reflection  of  these  poles  across  the  jw  axis.  Therefore,  the  original  poles 
can  be  found  by  choosing  only  those  left-half  plane  poles  from  the  set  of  all 
roots  of  (C15) . 

Note  that  up  to  this  point  in  the  derivation  only,  the  only  assumption 
used  is  that  the  function  is  causal.  Therefore,  the  poles  of  a causal  func- 
tion can  be  found  from  the  knowledge  of  magnitude-only  data. 

Once  the  coefficients  a are  known,  the  equation 

a 


N-l 

V aa(ju)“  = 0 (C16) 

a=0 


can  be  solved  for  the  unknown  roots  which  gives  the  zeros.  For  a positive- 
real  responses,  the  zeros  of  the  response  are  in  the  left-half  plane.  It  is 
obvious  from  (C8) , the  roots  of  (C,16)  will  contain  the  original  system  zeros 
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plus  a reflection  of  these  zeros  across  the  jo)  axis.  The  original  zeros  can 
be  found  by  choosing  only  those  left-half  plane  zeros  from  the  set  of  all 
roots  of  (C16) . 

Wien  the  poles  and  zeros  of  the  function  are  known,  the  residues  of  the 
function  remain  to  be  found.  To  obtain  an  expression  for  the  residues,  begin 
by  equating  Eq.  (Cl)  and  ( C 3)  to  get 


N-l 


N 

y -A- 

a=l  (s  " sa) 


n <■-  za) 


N 

n 

a=l 


(C17) 


(s  - sa} 


By  multiplying  both  sides  of  the  above  equation  by  (s  - s ) and  letting  s = s , 

q 

we  can  express  the  residues  as 


R 

q 


N-l 


n <8  - a ) 


N 

n;  <■,  - 
a=l  M 


(CIS) 


where  the  supscript  in  the  product  means  the  deletion  of  this  term  when  a = q. 

Note  from  (C6)  that  d,  = - \l a...  , and  the  term  a,T  has  been  found  in  the 

i y N-i  N-i 

inversion  of  (Cll)  . 

Before  the  numerical  results  of  this  method  are  discussed,  a possible 

limitation  due  to  the  algebraic  structure  of  the  method  should  be  presented. 

2 

In  a previous  report  where  a method  was  developed  to  extract  poles  and  resi- 
dues from  the  magnitude  and  phase  of  the  frequency  response  function,  it  was 
shown  that  the  defining  equation  used  to  find  the  coefficients  were  in  powers 

of  n.  The  similar  equation  for  the  magnitude-only  method  (presented  in  (CIO)) 

2 

is  in  powers  of  2n.  The  previous  method  was  shown  to  work  for  only  20-pole 
pairs  due  to  the  large  dynamic  range  of  the  scaled  system  of  linear  equations. 
Since  the  defining  Eq . (CIO)  is  in  power  of  2n  instead  of  n,  the  dynamic  range 
of  the  magnitude-only  method  should  be  twice  as  large  as  the  dynamic  range  for 
the  method  using  the  magnitude  and  phase  data.  Therefore,  the  magnitude-only 
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program  developed  In  this  appendix  should  be  expected  to  extract  the  poles  for 
response  that  have  up  to  only  10-pole  pairs. 

A program  based  on  the  above  was  written  and  tested.  The  linear 
Eq.  (Cll)  is  solved  for  the  coefficients  following  which  the  poles  are  found 
by  solving  for  the  roots  of  (C15),  and  the  zeros  are  found  by  solving  for  the 
roots  of  (C16) . The  residues  are  then  found  by  using  Eq . (C18) . Some  results 
obtained  from  this  numerical  program  will  not  be  discussed. 

The  program  was  first  tested  by  using  a 2-pole-pair,  positive-real 
function.  The  results  for  this  2-pole-pair  test  are  shown  in  Fig.  C-la.  This 
figure  shows  the  given  and  extracted  spectrum  and  the  isometric  plot  for  the 
given  and  extracted  pole-sets.  Also  shown  in  this  figure  are  the  tabulations 
of  the  given  and  extracted  poles  and  residues  sets.  It  is  obvious  from  this 
figure  that  the  technique  does  work. 

The  case  when  a 10-pole-pair,  not-positive-real  function  was  used  as 
input  data  is  shown  in  Fig.  C-lb.  From  this  figure  one  observes  that  the 
correct  poles  are  found,  but  the  correct  residues  were  not  extracted.  Since 
the  response  function  was  causal,  the  correct  poles  are  expected  to  be  found. 
The  residues  found  from  this  program  are  the  residues  that  fit  the  poles  and 
coefficient  of  (Cll)  in  a positive-real  sense.  This  test  shown  that  the 
method  can  be  used  to  extract  the  correct  poles  from  a causal  function  that  is 
not  positive  real. 

The  program  was  then  used  to  extract  poles  and  residues  from  the  numer- 
ical generated  data  for  the  driving-point  impedance  of  a 50-meter  wire  driven 
at  the  center.  The  procedure  failed  to  extract  the  correct  poles  and  residues 
from  the  given  data.  It  found  some  correct  poles;  but  it  also  found  many 
incorrect  poles.  A preliminary  investigation  into  this  extraction  procedure 
seem  to  indicate  that  the  difficulty  lay  in  the  root-finding  routine.  The 
present  root-finding  routine  is  of  the  type  known  as  Muller's  method.  This 
method  searches  in  the  complex  plane  for  a root  in  a rather  random  fashion, 
then  deflates  the  original  function  by  the  newly  found  root.  The  newly 
deflated  function  is  then  searched  for  a new  root.  It  is  suprising  that  this 

method  of  root  extraction  does  not  work,  since  it  has  been  used  with  great 

1.2 

success  in  the  previous  reports.  In  these  previous  tests,  the  function 
used  had  roots  only  in  the  left-half  plane,  but  now  the  function  has  roots  in 
both  half  planes.  This  abundance  of  roots  seems  to  be  confusing  the 
numerical-search  procedure  used  to  find  the  poles  and  zeros. 
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A new  numerical  root-finding  routine  was  obtained  that  finds  roots  by 

12 

method  developed  by  Jenkins  and  Traub.  This  method  searches  in  the  complex 

plane  for  the  root  with  the  smallest  modulus,  then  deflates  the  original 

function  by  the  newly  found  root.  The  newly  deflated  function  is  searched 

for  its  root  having  the  smallest  modulus.  The  procedure  used  to  search  for 

the  roots  in  this  method  is  different  from  all  previous  methods  and  mav  be 
13 

superior.  The  Jenkins-Traub  root-finding  routine  was  in  the  process  of 
being  substituted  for  the  Muller  routine  when  the  project  terminated. 

This  appendix  has  presented  a method  for  extracting  the  poles  from  mag- 
nitude of  causal  responses  and  the  poles  and  residue  from  magnitude  of  a posi- 
tion real  function.  This  method  should  be  valid  for  extracting  up  to  10  pole- 
pairs.  The  method  was  tested  and  shown  to  work  for  given  data,  but  failed  to 
work  on  numerically  generated  EM  data.  Tbe  latter  problem  apparently  stems 
from  the  method  by  which  the  roots  are  extracted  from  the  polynomial. 

In  the  future  the  new  Jenkins-Traub  root-finding  routine  should  be  sub- 
stituted for  the  routine  using  Muller's  method  to  find  the  roots. 


LRS/af / jn/vt/ j f 
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ANAI  YI  |f  FUNCT  ION  MAO 
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Fig.  C-la.  This  series  shows  the  results  from  testing  the  magnitude-only 
frequency-domain  Prony.  This  specified  transfer  function  along  the  ju 
axis  and  pole  location  in  the  o-ju  plane  (with  the  residue  value  indicated 
by  the  vertical  lines)  are  included.  The  data  shown  are  for  a 2-pole  pair, 
causal,  positive-real  function. 
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ANAL n ! < TUNC  T ION  MAO 


Spectrum  of  given  data 
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Fig.  C-lb.  Same  as  Fig.  C- la,  except  the  data  shown  here  are  for  a 10-pole 
pair,  causal,  not  positive-real  function. 
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